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I. INTRODUCTION 


Deuterium has been used in the past as a target in several neutrino quasielastic scattering 
experiments for average energies in the range 0.5 - 27 GeV, where the hnal-state particles 
in the reaction were detected in bubble chambers [T]. For reactions below pion production 
threshold, a neutrino interacts with a neutron bound in deuteron, which turns into a proton 
and a negatively-charged lepton is produced. The two resulting protons are no longer bound 
and one or both of them can be detected in coincidence with the charged lepton. 

From a theoretical point of view, the advantage of using deuterium as a target with coin¬ 
cidence detection lies in the fact that the kinematics of the remaining nucleon, and of all of 
the particles taking part in the process, are fully determined by energy-momentum conser¬ 
vation. This includes the energy of the incoming neutrino, which usually has a broad energy 
spectrum, an issue for neutrino oscillation experiments when nuclei other than deuterium 
are used. In the case of nuclei in general, the nucleus is in general excited to conhgurations 
at high missing energy (see [2] for a discussion of semi-inclusive CCz/ reactions including 
the dehnition of missing energy) and, even when both a charged lepton and a proton are 
detected in coincidence, one cannot reconstruct the incident neutrino energy. However, the 
deuteron does not share this problem and all kinematic variables can be reconstructed from 
measurements of a subset of hnal-state particles. The details of such procedures are dis¬ 
cussed later. We note that, insofar as a cut can be made to separate events where pions 
are produced from those where they cannot, i.e., where the “no-pion” cross section can 
be isolated (the kinematics for making such a cut are discussed below), what we continue 
to call the semi-inclusive cross section is actually an exclusive cross section, meaning that 
the energies and momenta of all particles are determined by measuring only the subsets 
summarized above. 

In addition to determining the neutrino energy using only the kinematics of the reaction 
being studied, namely, -|- ^H—> + p + p, the cross section for this reaction can be used 
to determine the neutrino hux. This, of course, requires that one knows that cross section, 
and indeed, especially under favorable conditions, this is the case as discussed later. Such is 
not the case, however, for complex nuclei such as carbon or oxygen where considerable effort 
has gone into evaluating the level of theoretical uncertainty in modeling neutrino reactions 
in those cases. 
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Aside from providing a means to determine the energy and flux of the incident neutrino in 
such semi-inclusive CCi^ reactions, neutrino disintegration of the deuteron has the potential 
to yield valuable new information on the nucleonic content in the problem. Specifically, 
given that the vector EM form factors of the nucleon are relatively well-determined from 
electron scattering on protons and light nuclei, the prime candidate for such an approach 
is the isovector, axial-vector form factor of the nucleon, Once the uncertainties in the 

modeling of the deuteron and pp systems have been evaluated, i.e., the extent to which 
different ground-state wave functions yield different results and to which the treatment of 
the pp hnal state gives different answers, the reaction can be used to determine Both 
aspects of the problem are discussed below. 


The paper is organized as follows: In Sect. [IT] the basic kinematics of the process is 
described, and the formalism of the weak responses and the neutrino-deuteron cross section 


is given in Sect. |III[ In Sect. the deuteron structure model is summarized. In Sect. 
[V| results are presented and discussed for inclusive and semi-inclusive neutrino-deuteron 


scattering for several choices of kinematics. Finally, in Sect. VI our conclusions are given. 


II. KINEMATICS 


Let us begin from a “theoretical” point of view and presume that the energy of the incident 
neutrino is known, having mass m with three-momentum k and total energy e = + m^, 

contained in the four-momentum = {e, k). Later we will work backwards, assuming only 
hnal-state momenta are known and show how the neutrino energy may be reconstructed. 
The four-momentum corresponding to the outgoing charged lepton, with mass m', is K''^ = 
(V,k'), and the four-momentum transfer is = (ti;,q), with —Q^ = |Q^| = > 0 

(spacelike). The three-momentum transfer q is assumed to be along the 3-axis so that the 
lepton momenta define the 13-plane (see Fig. [^, and the angle between them is the scattering 
angle 6. With these dehnitions the components of the above mentioned four-momenta can 
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FIG. 1. Kinematics for coincidence neutrino-deuteron reactions. 


be written as 


= e K'^ = e' = e - e' = u 

= -kk' sinO K'^ = -kk' sinO = 0 

= 0 K'^ = 0 = 0 

= ^k {k — k'cos 6 ) K'^ =—^k'{k'— kcos 6 ) = \^k‘^ + k'"^ — 2 kk' cos 6 = q . 

( 1 ) 

Concerning the hadronic part of the process, the deuteron at rest carries four-momentum 
= (Mrf, 0,0,0), where Md =1875.61 MeV is its mass. The protons resulting from the 
CCz/ disintegration of the target have four-momenta Pf = (Pi,pi ) and = (P 2 ,P 2 ) 
respectively. Let d[ and 62 be, respectively, the angles of the proton momenta with respect 
to the incident neutrino momentum and and 02 their angles with respect to the leptonic 
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plane; the components of the proton fonr-momenta can then be written as 


P? = E, 


Pi = P2 [— sin 6*2 < 
and the 




1 

Pi = Pi [sin 9[ cos (j)'i cos 9q + cos 9[ sin 9q] 

= Pi sin 9 'i sin 0 '^ 

' 16*g + cos 6*'j^ cos 6*q 

momentnm transfer q. 


Pi 

Pf — Pi [~ ®i^ ^1 CO® i^'i ^i'^ 
where 9 q is the angle between the incident momentnm k 


Pi — P2 sin 6*2 sin i 

6*2 cos (f)2 sin 9q + cos 6*2 cos 9q] , 

( 2 ) 


so 


that 

sin 6*0 = -fc'sin6* and cos6*„ = -(k — k' cos9) . ( 3 ) 

q q 

The primed angles 9[, 6*2, and dehned with respect to the incident nentrino mo¬ 

mentnm, are related to the nnprimed ones dehned with respect to the momentnm transfer, 
as were used in [2], through 


sin 6*' cos 0' = sin 6*^ cos 0j cos 9 q — cos 6*^ sin 9 q 

sin 9[ sin 0 ' = sin 6*j sin (j)i ( 4 ) 

cos 9[ = sin 6*j cos 0 * sin 9q + cos 6*j cos 9q 


with 02 = 01 + vr. Let us next turn to the reverse situation, the “experimental” point of 
view where only hnal-state particles are presumed to be detected and where the goal is to 
determine the incident neutrino energy. 


A. Neutrino energy determination 

On the one hand, energy and momentum conservation in the hadronic vertex implies 
u: = El + E2 — Md and q = Pi + P2, or g = pi cos6*1 ± y^pf^-pfsh?^ for its magnitude. 
On the other hand, from the lepton vertex one has, as shown in Eqs. Q, cu = £ — e' and 
q = \/k‘^ + — 2kk' cos 9. The transfer variables oj and q must be equal in both vertices, 

and positive. Hence, the energy-momentum conservation conditions yield two relationships 
between the undetected proton energy ^2, or corresponding momentum p2, and the incident 
neutrino energy e, or corresponding momentum k] the only possible values of E2 and £ (or 
P2 and k) for the process are those that fulhll both conditions simultaneously. It should be 
noted that this procedure cannot be followed with a general nuclear target, since the rela¬ 
tionship between E2 and p2 is not hxed. Indeed, although on-shellness always holds, namely 
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FIG. 2. (Color online) Left: An example of using energy conservation (dashed curve) and momen¬ 
tum conservation (solid curve), to determine the incident neutrino momentum k as a, function of 
the undetected proton momentum p 2 for a deuteron disintegration with emitted muon momentum 
k' = 1 GeV, scattering angle 9 = 60°, emitted proton momentum pi = 0.5 GeV, and proton emis¬ 
sion angle 9i = 20°. The intersection of the curves gives the only possible combination of incoming 
neutrino and undetected proton momenta for this process (1.45 GeV and 0.83 GeV, respectively). 
Right: Using this procedure, and for the same kinematic conditions for 6, pi and 0i, resulting 
incident neutrino momentum k (solid curve) and undetected proton momentum p 2 (dotted curve) 
as a function of the emitted muon momentum k'. 

E\_^ = p\_i + the rest mass of the residual system Wa-i is in general unknown (see 

0); however, in the specihc case of the deuteron target, Wa-i = fnp without ambiguity. 
Fig. a shows a graphical example of how the constraints work. In the following subsection 
we develop this possibility of neutrino energy determination in deuteron scattering analyt¬ 
ically by defining and studying specific kinematic scenarios, as well as the pion production 
threshold and other kinematic constraints. 

B. Kinematic scenarios 

We consider two scenarios, the first (‘scenario A’) where the two protons in the final 
state are presumed to be detected, but not the muon, and the second (‘scenario B’) where 
the final-state charged lepton and one proton (called proton number 1) are presumed to be 
detected, but not the other proton (proton number 2 ). Using the nomenclature given above. 
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we define the following fonr-momenta: 

P':, = (EA,p^) = P^ + P!^ ( 5 ) 

P^b = {Eb,Pb) = K"‘ + P; , ( 6 ) 

which yield the following relationships: 

Ea = El + E 2 = Md + CO ( 7 ) 

Pa=Pi+P 2 = <1 (8) 

Eb = s' -\- El ( 9 ) 

PB = k'+Pi. ( 10 ) 


Note that in scenario A this means that q and u: are immediately known. The angles between 
the incoming nentrino direction and the three-vectors and Pb are known and denoted 
9 a and 9 b, respectively. Referring to Fig. [^we see that 9 a = 9 q. Let us now develop the 
two scenarios one at a time. 


Scenario A: In this case one knows q and 00 through q = Pa and 00 = Ea — M^. Noting 
that k' = k — q, one has 

= k‘^ + — 2 kq cos 9 q (11) 

and together with s' = e — u, and using the fact that — k''^ = one obtains 


k q cos 9 q — 00 e = Xa 


( 12 ) 


with 

Xa = ^ - 1 - m'^ — > 0 . 

This immediately leads to values of the incoming neutrino momentum and energy: 


1 

Qa 

geos 9 qXA -l- 

! X\ + m?aA 

1 

Qa 

ojXa + q cos 9 q 

j X\ + m?aA 


where 


qa = {q cos 9 q)‘^ — > 0 . 


( 13 ) 


( 14 ) 

( 15 ) 


( 16 ) 
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One can show that O<6'q<0°<7r/2, where 0 ° = arccos(u;/g). Knowing the neutrino 
beam energy e then yields the charged lepton energy s' = e — u and thus its momentum 
k’ = together with the angle 6 through 



k' cos 9 = k — q cos 6q 

k' sin 6 = q sin 6q . 

(17) 

(18) 

That is, one now has all of the kinematic variables in scenario A. 


Scenario B: In this situation Eb and are presumed to be known, and the unknown 

proton variables can be expressed as E 2 = e — {Eb — Md) and P 2 = k — Pb', from the latter 

one gets 


pI = k'^ + p% - 2 k pb cos 9b , 

(19) 

and then using the fact that E 2 — p\= one has 



k pb cos 9b - {Eb - Md) e = Xb 

(20) 

with 

Xb 

from which one obtains 

= ^[pb- {Eb - Md)"^ + >0 

(21) 

k= — 

Ub 

Pb cos 9bXb + {Eb - Md)\jx% + m?aB 

(22) 

1 

Ub 

{Eb - Md) Xb + Pb cos 9B\jx% + m?aB 

(23) 

where 

aB = {pb cos 9b)‘^ - {Eb - Md)^ • 

(24) 


In this case, one now has q and u using the usual relationships in Eqs. ([^ and can then 
hnd E2 and Pg using Eqs. ([^ and (|^, yielding all of the kinematic variables in scenario B. 

We now study the incident neutrino energy or momentum threshold to produce the 
lightest possible extra particle in the scattering process, namely a neutral pion, vr^. One 
assumes the knowledge of the three emitted particles variables, and then dehnes 

= [EcPc) = K'f + P‘‘ + , 
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(25) 









so that Ec = e' + El + E2 and P(j = k' + + P2- One then has Et,o = Md + e — Ec and 

p^o = k — Pq, whence 

pIo = k'^ + Pc- 2 k Pc cos Oc , ( 26 ) 

and using the on-shell relation for an emitted neutral pion, E'^o — one obtains 

k Pc cos 9 c - {Ec - Md) e = Xc , ( 27 ) 


with 

= \ [Pc - {Ec - Mdf + mlo - m^] . ( 28 ) 

From these expressions one obtains the following pion production threshold values of the 
incident neutrino momentum and energy: 


k-j^o ih 


th. 


1 

ac 

1 

ac 


Pc cos Oc Vc ±(Ec - Md) Xc^mlac 
{Ec - Md) Xc ±pccosOc^Xc + m^ac , 


(29) 

(30) 


where 


ac = (pcaosOc)'^ - {Ec - Md)^ ■ 


(31) 


For given energies of the emitted particles e', Ei, E2, neutral pion production is ruled out 
in the scattering process if the incident neutrino beam energy fulhlls e < (and equiv¬ 

alently for the momenta). 

A potential strategy for measurements of the desired kinematic variables might be the fol¬ 
lowing: one might assume that the three particles in the hnal state, say, a muon and the 
two protons, are all measured with adequate precision hrst to eliminate the possibility of 
7 r° production, or, more generally, to isolate processes where a pion is produced from those 
where it is not, since the former is interesting in its own right. Given this hrst cut, one can 
then safely proceed to analyze the reaction 


p, +P + P 


as above in either of the scenarios, whichever proves to be the more favorable from an 
experimental point of view. 

Finally in this section, for completeness it is also useful to make contact with the general 
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developments of semi-inclusive CCz/ reactions presented in [2]. The variables introduced 
there translate into the present ones in the following way: 


O Mrf ; ruAT O rup ; Wa-i ^ nip ; Pjv ^ Pi : P^ -P 2 ■ (32) 

The momentum of what is labeled particle 2 above is usually called the missing momentum: 
Pm = P 2 = ~P- For simplicity below we will show results as functions of p = |p|. The 
quantity in that paper called T in previous scaling analyses (see ma) is zero, since the 
daughter system in the general case is simply the second proton in the present study. This 
leads to the following for the scaling variable y and the quantity V: 



(33) 

(34) 


Y = q + y, 


where the invariant mass in the final state is given by 



(35) 


and its threshold value by Wt = 2mp, so that W > Wt- A useful relationship that emerges 

is 



(36) 


and an important constraint in the reaction is the following: 


\y\<p<Y . 


(37) 


Furthermore the cosine of the angle 9pq in [2] is given by 



(38) 


III. WEAK RESPONSES AND CROSS SECTION 

Using the Feynman rules and integrating over the undetected nucleon three-momentum, 
the cross section of the process can be written as 



where is the matrix element squared (see below), uq = (^ + ~ ^Fe differentials 

can be expressed as = k'‘^dk'dVLk' = k's'de'dVLk' and d^pi = p^dpidflp-^ = piEidEidVLp^. 
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From this expression one can construct the differential cross section with respect to any set 
of variables, taking into account that the integration of the delta function 6{f{x ,...)) with 
respect to the variable x forces the variables of / to fulfill the condition / = 0 and introduces 
an extra factor \df{x, ...)/dx\~^. 

As an example, if the energy distribution of the incoming neutrino beam, P(e), is known, 
the differential cross section averaged over incident energies is given by 
f day 


de' dVtk' dEi dVt 




'Pi 


2(27r)' 


Pi ^0 p. . I p 


-1 




with 


Fp = 1 — 


p\ COS 01 

9 


) (1-^) 


(40) 


(41) 


As another example, by a further integration of Eq. (39) over the undetected nucleon 
energy one gets the following differential cross section: 


day 


dk' dVty dVtp^ 


cos^ 9c rrip k'"^ vq 
k e' 


2 

_p 


1-1 


with 


F = 1 + 


2(27r)5 
Lu Pi — q El cos 9i 




MdPi 


which becomes in extreme relativistic limit 
da X [erl] 


cos^ 6c ml X 1 


Ti 


(42) 


(43) 


(44) 


de' dVty dVtp^ 

The energy and momenta in these expressions fulfill conservation laws and thus the unde¬ 
tected energies e and E 2 (or momenta k and P 2 ) take only specific values that can be deduced 
as described in the previous section. 

The matrix element squared is given by the contraction of the leptonic and hadronic 
tensors, which can be written in terms of products of generalized Rosenbluth factors V and 
hadronic responses w for charge, longitudinal and transverse projections: 


+ Rcl + Vll + wtt) + 

Vt + Vtt {zOtt + zv^t) + ^TC {ziZtc + ^Tc) + 

+Vtl {zvtl + zv^l) + X [Vr'zv'^f^ -1- Vtc'Zv^c' + ^tl'Zv^i,] , 


(45) 


where y = -|-1 for neutrino scattering and x = — 1 for antineutrino scattering. The gen¬ 
eralized Rosenbluth factors contributing to this matrix element are given in [2] for general 
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lepton masses using the quantities 5 = mj and b' = m'/ a/IQ^I and for general val¬ 
ues of the vector and axial-vector coupling constants ay and a^- Here we reproduce those 
expressions particularized to the Standard Model tree-level values of the coupling constants 
{ay = —cla = 1) and to (anti)neutrino scattering with m = 0: 


Hcc 

VcL 

Vll 

Vt 

Vrprp 

Vtc 

Vtl 

Vj'i 

Vtc 

Vtu 


1 — 5'^ tan^ 9/2 
—v - - 6'^ taiV 9/2 


P 


9 / 

7 

2u 


1 + 

— 

\ 

V 

P' 

1 

/ 

V 

2^ + 

V- 

P' 


s:'2 


/2 \ r /2 


r/4 


- 2 '’ 


1 + - + 


■/2 


6'^ tan" 0/2 


—- tan 9/2 \ /— Vtt 
P VP 

h> -f pp' (5'") Vtc 


- \- p S'^ tan" 9/2 

P ' 


tan 0/2 - 


P L 


' + 7 + 2 '’" 


/2 


tan" 0/2 


1/2 


-vVtc' ■ 


(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 


The following definitions have been used in the expressions above: 

(se) 

The deuteron responses w, on the other hand, can be constructed as functions of the three 
independent four-momenta of the hadronic part of the process, namely Pj/, and P^, in 17 
different VV, AA and VA combinations, each associated to an invariant function of the four 
dynamical invariants constructed with the same set of four-momenta [2]. The computation 
of these deuteron responses will be described in the next section for a sophisticated model 
of the nucleon structure. 


IV. DEUTERON STRUCTURE 

The calculation of matrix elements for electro-weak breakup of the deuteron requires 
that some care be exercised in the construction of a consistent model of the reaction in 


12 









order that basic symmetries, such as electromagnetic current conservation, be maintained. 
This is a problem that has been discussed by many authors in the context of deuteron 
electrodisintegration. A detailed discussion in the context of the Bethe-Salpeter equation 
[6] and the related spectator equation [TJ ^ can be found in |9l HO]. In general, for all 
approaches a nucleon-nucleon [NN) interaction is parameterized in terms of a potential or 
interaction kernel which can be iterated in the appropriate equation, such as the Schrodinger 
equation or Bethe-Salpeter equation, to produce the deuteron bound-state wave function and 
scattering amplitudes which is then fit to data for laboratory kinetic energies up to Tiab ~ 
350 MeV, which is slightly above pion production threshold. Since the NN interaction 
generally involves the exchange of electro-weak charges, it is necessary in constructing the 
model to include two-body currents that contain coupling to the exchanged particles that 
carry the charges. This program generally relies on the ability to relate the NN interaction 
to some set of effective meson exchanges and has been carried out rigorously, for example, in 
the context of chiral effective held theory m and the covariant spectator equation [T^HTT] . 
However, attempts to extend this approach to invariant masses well above pion threshold 
have been relatively unsuccessful and at large four-momentum transfers and invariant masses 
it is necessary to construct models which are not fully consistent as a result of the inability 
to produce interaction kernels that can successfully describe NN scattering in this region. 

The calculations presented in this paper use a model related to the Bethe-Salpeter equa¬ 
tion, which is designed to be used at large four-momentum transfers and invariant masses 
and has been used previously to describe deuteron electrodisintegration [T5tiT9] . A brief 
outline of the model containing the modihcations required for calculation of the CCu reac¬ 
tion will be presented here. This model is constructed primarily for use at large where 
relativistic effects are important and where there are open meson-production channels that 
must appear as inelasticities in the scattering matrix producing the hnal-state interaction. 
However, the plane-wave contributions can be used at all values of Q^. At this stage no 
two-body meson-exchange currents (MEC) are included in the calculations, although these 
will be included in the future, meaning that the calculations shown here are in the impulse 
approximation (lA). As a rough estimate of what this approximation might entail, we note 
that in [20] the scaling behavior of EM (vector) MEC was found to be such that the results of 
the ratio between MEC and lA goes as k\, where kA is some characteristic three-momentum 
for a given nucleus (roughly the Fermi momentum kp). This ratio is typically found to be 
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FIG. 3. (Color online) Feynman diagrams representing charge-changing neutrino reactions in the 
impulse approximation. 

about 15-20% at the maxima of the quasielastic and MEG contributions for nuclei such as 
where = 228 MeV. However, the characteristic momentum for is about 55 MeV, 
which yields a rough estimate for this ratio in the present situation of ~5-7xl0“'^, namely, 
contributions from MEG should be very small. As a result, this calculation is limited to the 
lA. 

The Feynman diagrams representing the impulse approximation are shown in Fig. 

Fig. 1^ represents the direct plane-wave contribution. Fig. represents the corresponding 
exchange contribution and Fig. represents the contribution from hnal-state interactions. 
Use of Fig. alone is called the plane-wave impulse approximation (PWIA), while adding 
the exchange diagram Fig. gives the plane-wave Born approximation (PWBA) [211 
122] • We refer to the inclusion of all three diagrams in Fig. as the distorted-wave Born 
approximation (DWBA). In the next section results are shown for all three assumptions to 
ascertain the impact of antisymmetrization and hnal-state interactions on the predictions 
made in this study. 

The matrix element corresponding to is given by 

(pisi; P 2 S 2 I J^{Q) |pdAd)„ = -h(pi, si)P^c(g)Go(Pd - P2)r^^(h^2, Pd)h^(p2, S 2 ), (57) 

where the target deuteron has four-momentum and spin A^, the two hnal-state protons 
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have four-momentum and spin (Pi,si) and {P 2 ,S 2 ). The single-nucleon propagator is 


G,(P)= 


— P'^ — if] 

and the weak charged-current operator has vector and axial-vector contributions 


( 58 ) 


^^cciQ) = KiQ) - 


KiQ) = 


n{Q) = 


Ga(QV + g'p(q 


2m jq 
2^ 


2m N 


75: 


(59) 

(60) 
(61) 


where the vector-isovector form factors used here are from the GKex05 parameterization (see 
[22] and references therein) of the nucleon electromagnetic form factors and the axial-vector 
and induced pseudoscalar form factors are taken to be 


Ga{Q") 


1.2695 



(62) 


and 

Gp(Qt = 

f _L _21 ] 

y 185.05 ^ 4m2 J 

Unless stated otherwise we use Ma = 1.03 GeV. 

For this contribution, the one leg of the deuteron is on-shell. The deuteron vertex function 
with nucleon 2 on-shell can be written as 


rA,(P2, Pd) = 9l{Pl P2 ■ Pd)7 ■ ^xAP) + 92{PI P2 • Pd) 




g,{Pl P2 • P,)7 ■ ^x^Pd) + 94{PI P2 • Pd) 


P-^xAPd) 


rriN 


rriN 

7 ■ Pi + m 

rriN 


C, 


(64) 


where Pi = P^ — P 2 , P = |(Pi — P 2 ) = ^ ~ P 21 C is the charge-conjugation matrix and 
is the deuteron polarization four-vector. The invariant functions 9 i are given by 


9i{PlP2-Pd) = 


92{PlP2-Pd) = 


9,{PlP2-Pd) = 


2Pk — Md 
■y/S^ 

2Pk — Md 

■y/s^ 




3 mN 


mN , I , mNi2E^ + mN) , . , /3 tuat 

-u{k) H-- w{k) + \ - - Vt{K) 


E^ + mN 


V2, 




2 K 


3 mNEf, 
16n K 


Vt{P) 


(65) 

( 66 ) 

(67) 


15 

























■ Pd) — 


m 


N 


y/MMd 


(2-E/k — Md) 


Ek + ruN 


-U[K] 


+ 2mN , A , \PiMd . - 
-w{k) H- Vs{k) 


V2, 




K 


( 68 ) 


where 


K = 


\Pd-P2)^ 

PI 


PI 


is the magnitude of the neutron three-momentum in the deuteron rest frame and 


= \! k'^ + mj,. 


(69) 


(70) 


The functions m(k), w{hi), Vs{E) and Vt{E) are the s-wave, d-wave, singlet p-wave and triple 
p-wave radial wave functions of the deuteron in momentum space. For convenience, the 
spectator deuteron wave function can be defined as 


i’\d,S2{P‘2-, Pd) — Go{Pd — P2)TT{P2, Pd)u^{P2, S 2 ) 


(71) 


We choose to normalize this wave function such that in the deuteron rest frame 


E f Pdh‘‘ix..,iP2,Pi) = I, 


S2 


(27r)3 E.^ 


(72) 


which is correct only in the absence of energy-dependent kernels |9]. This corresponds to a 
normalization of the radial wave functions 


dKK^ 


3 ['u“(k) + y/{K) + v^(k) + !I3(k)] = 1 . 


Jo (2!r) 

For the calculations in this paper we use the WJC 2 wave functions 

The current matrix element corresponding to Fig. can then be written as 

{pisi;p2S2l J^iQ) |pdAi)„ = -ti(pi.si)rg(,((())ji>j3,„(P2.7)i). 


(73) 


(74) 


The current matrix element corresponding to Fig. ^ is related to that of Fig. ^ by a 
crossing of the two hnal-state protons. So, 


(pisppsSsI J^iQ) IPdAd);, = (P2S2;PiSi| J^{Q) IPdAd)^ • 


(75) 


The contribution from final-state interactions represented by Fig. requires the intro¬ 
duction of a pp scattering amplitude M and an integration for the loop four-momentum 
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k 2 , which involves both the deuteron vertex function and the pp scattering amplitude. In 
this case, both of the protons are in general off-shell. However, an examination of the con¬ 
tributions of the poles of the nucleon propagators shows that the contribution from the 
positive-energy on-shell pole of proton 2 dominates the calculation. As in the previous elec¬ 
trodisintegration calculations, we choose to put particle 2 on its positive-energy mass shell to 
simplify the calculation. Using this approach, the contribution of the hnal-state interaction 
to the current matrix element is given by 


(piSi;P2S2| J^{Q) \Pd>^d)c 


d?k 


2 m 


Ma(pi, Si)Ub{p2, S2)Mab-cd{Pl, P 2 ] K 2 ) 


J ( 2 vr) 3 Ufc, 

XGoceiPd + Q- K2)T^ccefiQ)GofaiPd - K 2 ) 

xA.\(k2)rI^^,(i^2,P), 


(76) 


where M is the pp scattering amplitude, 

A+(p) = “(P’ ^)“(P’ 

S 

is the positive energy projection operator and the Dirac indices for the various components 
are shown explicitly. 

At this point, only the incoming momentum for particle 1 is off-shell. For small values 
of the final-state invariant mass, this poses no problem, since the spectator equation could 
be used to construct a scattering matrix consistent with the deuteron bound state and a 
consistent current operator could also be constructed. However, for invariant mass well above 
pion threshold, there are no existing particle-exchange models of the kernel that reproduce 
the data. For this reason, it is necessary to use scattering matrices that have been £t to data, 
which limits the scattering amplitudes to the case where all legs of the scattering matrix 
are on-mass-shell. In na, the separation of the propagator for particle 1 into on-shell and 
off-shell contributions is described in some detail. All contributions can be calculated if a 
prescription is provided for taking the initial four-momentum of particle 1 off-mass-shell. 
This problem has been discussed in [23]. In the calculations presented here, we will use only 
the on-shell contribution which is well determined and is dominant. 

The completely on-mass-shell scattering amplitude can be parameterized in terms of hve 
Fermi invariants as 


Pkab\cd P S'('^5 t'}^ac^bd T ^)Tac ' bftd T .Fj'(s, i)^ac i^bii'')bd 
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( 78 ) 


where s and t are the usual Mandelstam variables and amplitudes Ti can be separated into 
isoscalar and isovector contributions using 



(79) 


The invariant amplitudes have been constructed in two ways. First, the helicity amplitudes 
obtained from SAID [23] - [?F] have been used to obtain the invariant functions as described in 


[T5] . In this case, the analysis of the pn amplitudes is limited to s < 5.97 GeV^ and the pp are 
limited to s < 9.16 GeV^. The second method is the use of a Regge model of NN scattering 
£t to cross sections and spin observables over the range 5.4 GeV^ < s < 4000 GeV^ [2811^ . 


With this brief summary of the dynamics used in the present study, let us now proceed to 
the presentation of some typical results both for inclusive and semi-inclusive GGz/ reactions 
with deuteron. 

V. RESULTS 

A. Inclusive cross sections 

We start by considering results for the inclusive cross section, i.e., where only the final- 
state charged lepton is presumed to be detected. We assume that only the no-pion cross 
section is allowed and that other channels are not involved. This, of course, is an issue for 
experimental studies and any inclusive measurement must involve all open channels. Our 
purpose in starting with modeling of the inclusive no-pion cross section is to ascertain where 
in q, u and scattering angle one should expect to find significant CCu strength before going 
on to study the semi-inclusive reaction. We begin by showing in Figs. and the inclusive 
no-pion CCu deuteron disintegration cross sections as a function of the momentum of the 
outgoing muon for three incident neutrino momenta {k = 0.5, 1 and 3 GeV) and for four 
small lepton scattering angles, 9 = 2.5°, 5°, 10°, and 15° (Fig. and four larger angles, 
6 = 45°, 90°, 135°, and 180° (Fig. [^. In each plot the corresponding momentum transfer at 
the peak (in MeV) is indicated, and results for PWIA, PWBA and DWIA calculations are 
shown. From these curves we extract the muon momentum at the peak of the inclusive no- 
pion cross section (/cq), as well as at one order of magnitude below the peak before and after 
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/cq ik'^ and respectively). As can be seen in these plots, the general trend is a decrease 
of the cross section when including exchange terms, he., from PWIA to PWBA, as well 
as when considering hnal-state interactions, he., from PWBA to DWBA. For momentum 
transfers around 500 MeV and larger the PWIA and PWBA results are indistinguishable; 
and they are also indistinguishable from the DWBA result for momentum transfers between 
500 and 1000 MeV, approximately. From the last we conclude that, at least for inclusive 
no-pion scattering when at the peak of the cross section and when the momentum transfer 
is reasonably large, the details of the modeling (he., as reflected in the three models used in 
the present work) are not too important. In contrast, for kinematics where the exchanged 
momentum is small, there are signihcant differences in going from model to model, and for 
quantitative comparisons with data some care should be exercised. 


B. Semi-inclusive cross sections 


For each combination of initial and hnal lepton momenta, k = 0.5, 1 and 3 GeV and 
k' = /cq, k'y (the latter dependent on k and 6*), we show in Figs. [^andj^CCz/ deuteron 
disintegration semi-inclusive cross sections as a function of the missing momentum (residual 
nucleon momentum) for two representative lepton scattering angles, 9 = 10° and 135°, and 
one representative nucleon emission angle (j) = 0°. We show PWIA, PWBA and DWBA 
results, which for 6 = 0° follow the same trend as for the inclusive no-pion cross sections, 
namely decreasing values in the order PWIA, PWBA, DWBA. For 6 = 135° the three 
calculations are very close. 

In Fig. 1^ we show semi-inclusive neutrino-deuteron cross sections for different nucleon 
emission angles 0 for representative neutrino momenta and scattering angle in DWBA at 
the quasielastic peak (k' = k'^). In order to clearly see the deviation of the cross sections due 
to the interference responses (the 0-dependent terms) that contribute to the cross section in 
DWBA, we show in Fig. [^for the same kinematics at the peak and in Fig. 10 at k' = k'^, the 
ratios of semi-inclusive neutrino-deuteron cross sections for different nucleon emission angles 
0 over the result for 0 = 0°. The deviations increase with increasing angle 0, increasing 
missing momentum p and decreasing incident momentum k, reaching up to a 40% at the 
quasielastic peak and up to a 60% at k' = k'^. Upon integrating to get the inclusive no-pion 
cross section, all of these 0-dependent contributions disappear; however, for semi-inclusive 
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FIG. 4. (Color online) CCu deuteron disintegration inclusive no-pion cross section as a function 
of the outgoing muon momentum k' for three incident neutrino momenta {k = 0.5, 1 and 3 GeV) 
and four lepton scattering angles {9 = 2.5°, 5°, 10°, and 15°). Dotted lines correspond to PWIA 
results, dashed lines to PWBA and solid lines to DWBA. Numbers in the top left corner of each 
plot correspond to the momentum transfer q in MeV at the peak of the cross sections. 
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FIG. 5. (Color online) As for Fig. but now for larger scattering angles {9 = 45°, 90°, 135°, and 
180°). 


measurements there are clearly significant effects to be seen when comparing in-plane with 
out-of-plane detection of the final-state protons. 


For completeness we have also studied the effect of initial-state interactions by using 


different deuteron wave functions. As an example we show in Fig. 11, for representative 
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FIG. 6. (Color online) CCi^ deuteron disintegration semi-inclusive cross section as a function of 
the missing momentum p for three incident neutrino momenta {k = 0.5, 1, and 3 GeV), and three 
outgoing muon momenta {k' = k'^, kg, and k'^ , see text). The lepton scattering angle is fixed to 
6 = 10° and nucleon emission angle to ()> = 0°. Results correspond to PWIA (dotted line), PWBA 
(dashed line) and DWBA (solid line). 


kinematics, the CCz/ deuteron disintegration semi-inclusive cross sections using different 
deuteron wave functions. For the kinematics shown, this set of wave functions yields results 
that differ by only a small amount. The largest deviation from the average is for the WJC 1 
deuteron vertex function which differs for the rest at p = 0 by less than 5%. The interaction 
kernel in this case contains a small admixture of pseudoscalar pion exchange which produces 
a hard short-range force that causes cross section strength to be transfered from small to 
large p resulting in an unusually small cross section at p = 0. 
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FIG. 7. (Color online) As for Fig. but now for lepton scattering angle Q = 135°. 


Effect of axial mass 


The above results have all been obtained using an axial mass Ma = 1.03 GeV in the 


parameterization of the axial form factor of Eq. (62). We have computed the same inclusive 


no-pion cross sections using an increased value of the axial mass, Ma = 1.3 GeV, and we 
have found that for small momentum transfers (g ^ 150 MeV), the results hardly change as 
expected, since then the Q^-dependence in the form factors is small and they are determined 
by their 0 limit. For larger momentum transfers, however, the Q^-dependence is now 

important and one sees that the increased axial mass gives noticeably larger inclusive no-pion 
cross sections. For momentum transfers between 150 MeV and 500 MeV, approximately, 
there is an overlap between the spread of results for different calculations (PWIA, PWBA 
and DWBA) obtained with the two values of the axial mass. Therefore, in that kinematic 
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FIG. 8. (Color online) Semi-inclusive neutrino-deuteron cross sections for different nucleon emission 
angles (p within DWBA for three neutrino momenta k = 0.5, 1, and 3 GeV, scattering angle 9 = 
45° and muon momentum k' = k'^ (at the quasielastic peak). 


region the effect of exchange and/or distortion can be confnsed with a modihed value of the 
axial mass, which calls for extra caution when interpreting experimental data. For instance, 
inclusive no-pion neutrino-deuteron measurements at incident momentum fc = 0.5 GeV and 
scattering angle 9 = 45° (corresponding to g ~ 360 MeV at the peak), compatible with an 
axial mass Ma = 1.03 GeV when interpreted through a PWIA calculation, would actually 
correspond to a larger axial mass, close to Ma = 1-3 GeV, if they were (more correctly) 


interpreted using a DWBA calculation. This comparison is shown in Fig. 12 (results for 
Ma = 1.03 GeV were also shown in the hrst plot of Fig. |^. 

As said above, the three types of modeling give very similar results at momentum transfers 
between 500 and 1000 MeV. The negligible spread seen there might lead one to naively 
think that those kinematic conditions are well-suited to determining the axial mass from 
experiment, as long as data from other regions are ignored in the analysis to avoid confusion. 
We show in Fig. [T^ examples of results in that region for Ma = 1-03 GeV and for Ma = 
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FIG. 9. (Color online) Ratios of the semi-inclusive neutrino-deuteron cross sections (as.-i.) for 
different nucleon emission angles (p over the result for </> = 0°, within DWBA for three neutrino 
momenta k = 0.5, 1, and 3 GeV, scattering angle 0 = 45° and muon momentum k' = fcg (at the 
quasielastic peak). 


1.3 GeV: at incident momentum k = 3 GeV and scattering angle 6 = 10°, corresponding to 
q ~ 530 MeV at the peak (left), and at k = 0.5 GeV and 6 = 135°, q ~ 700 MeV (right). 
However, even if the adjustment of the data to a given axial mass could be done in principle 
in that region, the corresponding momentum transfer is uncertain within a range related 
to the incident neutrino flux. As we showed above, only semi-inclusive measurements, not 
inclusive ones, allow one to determine the incoming neutrino momentum and therefore the 
momentum transfer of the process (within the experimental uncertainties). Semi-inclusive 
measurements and calculations are thus required to precisely determine the axial mass. We 


show in Fig. M, at the same kinematic conditions of the plots in Fig. IB and at the 
quasielastic peak {k' = /cq), the comparison of semi-inclusive cross sections for different axial 
masses. Glearly, measuring such semi-inclusive cross sections has the potential of yielding 
an excellent way to determine the axial-vector form factor. The uncertainties from the 
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FIG. 


10. (Color online) Same as in Fig. [^but for muon momentum k' = k'^. 


vector form factors from electron scattering are relatively small, and, as we have seen, the 
nncertainties from modeling the initial- and final-state two-body problem is arguably the 
smallest to be found in the entire periodic table. 


VI. CONCLUSIONS 

We have defined the kinematics of deuteron disintegration by neutrino scattering, + 
^H—)■ + p + p, and analyzed different detection scenarios in terms of both theoretically 
natural and experimentally practical kinematic variables. The deuteron as a target has the 
peculiarity that a semi-inclusive (coincidence) measurement of any two of the three particles 
in the hnal-state is actually exclusive, since the hnal state, consisting of a charged lepton and 
two protons (below the pion production threshold), has no other open channel, in contrast 
to the usual situation for complex nuclei where both missing momentum and missing energy 
dependences occur. One obvious conclusion is that semi-inclusive studies of this reaction 
hold promise for determining the incident neutrino energy using kinematics alone. 

In this study we have presented results for both inclusive no-pion and semi-inclusive 
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FIG. 11. (Color online) Semi-inclusive neutrino-deuteron cross sections using the deuteron wave 
functions AV18 [30], CD Bonn |3l|, N3L0500 [32|, Nijmegen 1, 2 and 3 [33|, and WJC 1 and 2 |T2] 
for fc = 3 GeV, 9 = 45° and k' = k'^. 


CCv cross sections using for the underlying dynamics a relativistic model of the deuteron 
and hnal-state pp structure that involves an approximation to the Bethe-Salpeter equation. 
We note that such an approach is completely covariant: the kinematics, the initial-state 
deuteron, the current operators and the hnal-state pp system are all fully relativistic, in 
contrast to much of the modeling of CCu reactions in general, where approximations are 
often employed (for a discussion of the role played by relativity for such reactions see [M]). 
We have provided inclusive no-pion and semi-inclusive CCz/ deuteron disintegration cross 
sections for a variety of kinematical conditions within PWIA, PWBA (he., introducing hnal- 
state hadronic exchange terms) and DWBA (he., using hnal-state hadronic interactions 
including exchange). The calculation of the two-body dynamics with the model used in this 
work, particularly when including exchange terms and hnal-state interactions (DWBA), is 
computationally demanding and therefore not well suited as direct input to neutrino event 
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FIG. 12. (Color online) For axial masses Ma = 1-03 GeV (dark curves) and Ma = 1-3 GeV (lighter 
curves), CCz^ deuteron disintegration inclusive no-pion cross section as a function of the outgoing 
muon momentum k' for incident neutrino momentum k = 0.5 GeV and scattering angle 8 = 45°. 
Dotted lines correspond to PWIA results, dashed lines to PWBA and solid lines to DWBA. 




FIG. 13. (Color online) Same as in Fig. 12 but for A; = 3 GeV, 9 
GeV, 6 = 135° (right). 


10° (left), and for k = 0.5 









FIG. 14. (Color online) For axial masses Ma = 1.03 GeV (dark curves) and Ma = 1.3 GeV 
(lighter curves), CCv deuteron disintegration semi-inclusive cross sections as a function of the 
missing momentum p for the same kinematic conditions of Fig. [T^ namely k = 3 GeV, 0 = 10° 
(left), and k = 0.5 GeV, 9 = 135° (right), and the muon momentum at the peak, k' = k^. Dotted 
lines correspond to PWIA results, dashed lines to PWBA and solid lines to DWBA. 

generators. However, we anticipate being able to provide simple parameterizations of the 
cross section in work that is currently in progress, and these will soon be made available for 
event simulations. 

At small and large momentum transfers {q ^ 500 and q ^ 1000 MeV) the approximation 
used for the interaction (PWIA, PWBA, DWBA) plays a significant role. Clearly anti- 
symmetrization (going from PWIA to PWBA) is seen to be important for such kinematics, 
with a somewhat smaller effect seen to arise from hnal-state interactions (going to the full 
DWBA). On the other hand, initial-state variations using state-of-the-art deuteron wave 
functions are observed to be minor. In contrast, for intermediate values of g, the three 
types of approximation lead to essential a universal result, showing very small sensitivity 
to both initial- and hnal-state physics. Thus, with an appropriate choice of kinematics, 
namely where the Q^-dependence of the axial-vector form factor plays a minor role (see 
below), results of semi-inclusive CCz/ disintegration of deuterium could not only provide the 
incident neutrino energy, as stated above, but could also yield the incident neutrino hux 
with relatively minor uncertainties from the modeling. 

Finally, in addition to the motivations for studying the -|- ^H—)■ + p + p reaction 
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given above, it can also serve as a way to improve our knowledge of the Q^-dependence of the 
isovector, axial-vector form factor. We have seen that the impact of changing the axial mass 
Ma used in a dipole parameterization of this form factor is signihcant in regions where the 
dependence on the modeling of the two-body problem is weak. The Myi-dependence increases 
for increasing momentum transfer, being negligible aX q ^ 150 MeV. We have shown that 
there is a range of intermediate momentum transfers (150 ^ q ^ 500 MeV) where the 
spreads of the curves for different approximations (PWIA, PWBA, DWBA) computed with 
axial masses ~ 1 and Ma ~ 1-3 GeV overlap, leading to a potential confusion in the 
extraction of the axial mass from experimental data unless the full model is employed. We 
close by noting that at some level inclusive measurements always suffer from their inability 
to be performed at a specific value of q, whereas the ability to determine all kinematic 
variables using semi-inclusive measurements, including the momentum transfer, provides a 
unique tool for studies of this type. 
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